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We show that continuous filling or wedge-wetting tran- 
sitions are possible in 3D wedge-geometries made from (an- 
gled) substrates exhibiting first-order wetting transitions and 
develop a comprehensive fluctuation theory yielding a com- 
plete classification of the critical behaviour. Our fluctuation 
theory is based on the derivation of a Ginzburg criterion for 
filling and also an exact transfer-matrix analysis of a novel 
effective Hamiltonian which we propose as a model for wedge 
fluctuation effects. The influence of interfacial fluctuations 
is shown to be very strong and, in particular, leads to a 
remarkable universal divergence of the interfacial roughness 
£± ~ (Tf — T) -1 ' 4 on approaching the filling temperature Tf, 
valid for all possible types of intermolecular forces. 
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There are two reasons why it is extremely difficult to 
observe interfacial fluctuation effects at continuous (crit- 
ical) wetting transitions in the laboratory Firstly, 
critical wetting is a rather rare phenomenon for which no 
examples are known for solid-liquid interfaces and only a 
limited number for fluid- fluid interfaces (^||. Secondly, 
the influence of interfacial fluctuations in three dimen- 
sions (d=3) is believed to be rather small jlj. For exam- 
ple, for systems with long-ranged forces, the divergence of 
the wetting layer thickness I on approaching the wetting 
temperature T w is mean- field-like, £~ (T w — T) _1 , and 
the only predicted effect of fluctuations is to induce an 
extremely weak divergence of the width (roughness) £j_ 
of the unbinding interface: £j_ ~ V — ln(T K , — T). Non- 
classical critical exponents and an appreciable interfacial 
width are onlypredicted for systems with strictly short- 
ranged forces B , but even here the size of the asymptotic 
critical regime is very small and beyond the reach of cur- 
rent experimental and simulation methods |3|||[| ■ 

The purpose of the present article is to show that these 
problems do not arise for continuous (critical) filling or 
wedge-wetting transitions |t]-|9| occurring for fluid ad- 
sorption in three-dimensional wedges. First, we show, 
contrary to previous statements in the literature ||] , that 
critical filling can occur in systems made from walls that 
exhibit first-order wetting transitions. Consequently, the 
observation of critical filling transitions is a realistic ex- 
perimental prospect. Second, we argue that interfacial 
fluctuations have a strong influence on the character of 
the filling transition and, in particular, the interfacial 
roughness of the unbinding interface, which is shown to 
diverge with a universal critical exponent. The fluctu- 
ation theory we develop is based on the derivation of a 



Ginzburg criterion for the self-consistency of mean-field 
(MF) theory and also an exact transfer matrix analysis of 
a novel interfacial Hamiltonian model for wedge wetting 
which we introduce to account for the highly anisotropic 
soft- mode fluctuations. This model leads to a complete 
classification of the critical behaviour in d = 3 and pre- 
dicts some remarkable fluctuation dominated phenomena 
which we believe may be tested in the laboratory. 

To begin, we recall the basic phenomenology of wedge- 
wetting and highlight the mechanism by which critical 
filling occurs in wedge geometries even for walls exhibit- 
ing first-order wetting transitions. Consider a wedge (in 
e? = 3) formed by the junction of two walls at angles ±a 
to the horizontal (see Fig. 1). Axes (x,y) are oriented 
across and along the wedge respectively. We suppose the 
wedge is in contact with a bulk vapour phase at tempera- 
ture T (less than the bulk critical value T c ) and chemical 
potential /x. Macroscopic arguments |?],|| dictate that at 
bulk coexistence, /i = /i sat (T), the wedge is completely 
filled by liquid for all temperatures T c > T > Tp where 
Tp is the filling temperature satisfying Q(Tp) — a. Here, 
0(T) is the temperature dependent contact angle of a 
liquid drop on a planar surface. Thus, filling occurs at 
a temperature lower than the wetting temperature T w 
and may be viewed as an interfacial unbinding transi- 
tion (of first- or second-order) in a system with broken 
translational invariance. We refer to any continuous fill- 
ing transition occurring as T — > Tp, /i — > Hsat(Tp) as 
critical filling. Also of interest is the complete filling 
transition which refers to the continuous divergence of 
the adsorption as \i — > [i sa t (T) for T c >T>Tp which is 
known to be characterised by geometry dependent crit- 
ical exponents |l(J. Here, we focus exclusively on criti- 
cal filling and, in particular, the critical singularities oc- 
curring as t = (Tp — T)/Tp — > + at bulk coexistence. 
The phase transition is associated with the divergence 
of four lengthscales (see Fig. 1) each characterised by a 
critical exponent: the mid-point (x = 0) height of the 
liquid- vapour interface Iq the mid-point interfa- 

cial roughness £j_ ~ t~ u± , the lateral extension of the 
filled region £ x ~ t~ v * and the correlation length of the 
interfacial fluctuations along the wedge £ y ~ t~ u f. So 
far, there has been no discussion of the values of these 
critical exponents for three dimensional systems beyond 
a simple MF calculation for Iq ||. On the other hand, 
transfer-matrix studies Q in d — 2 indicate that fluctu- 
ation effects are very strong at wedge-wetting and lead 
to universal critical exponents f3o = ux = u x = l. This is 
highly suggestive that fluctuation effects play an impor- 
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tant role in d = 3, relevant to experimental studies. 




FIG. 1. Schematic illustration of an interface configuration 
in the wedge geometry showing the relevant diverging length- 
scales at the filling transition. The planar adsorption l v and 
planar transverse correlation length fn remain finite at the 
transition. 

Previous MF analysis || have shown that a suitable 
starting point for the study of wedge-wetting in open 
wedges (small a) is the interfacial model 



H\£\ = 



dx dy 



W{i-a\x\) 



(1) 



where £{x, y) denotes the local height of the liquid-vapour 
interface relative to the horizontal, £ is the liquid-vapour 
interfacial tension and W(£) is the binding potential 
modelling the wetting properties of the wall. At MF 
level, this functional is simply minimised to yield an 
Euler-Lagrange equation for the y-independent equilib- 
rium height profile £(x), S£ = W'(£ — cv|a;|), where the 
dot and the prime denote differentiation with respect 
to x and £ respectively. This differential equation is 
solved subject to the boundary conditions £(0) = and 
£{x)— a\x\ — as >oo. Here, £„ denotes the MF pla- 
nar wetting film thickness (i.e., W '(l n ) = 0) and remains 
microscopic at the filling transition. Integrating once 
the equation yields a simple equation for the midpoint 
height, Ea 2 / 2 = W(£ ) - W(£ w ), which can be solved 
graphically Note that at bulk coexistence, Young's 
equation implies W(£„) = — £6 2 /2 (within the present 
small angle approximation) so that the present model 
immediately recovers the macroscopic result Q(Tp) = a. 
Depending on the form of W(£) (at Tp) the divergence of 
£q as T — >Tp is first-order or continuous. The condition 
for critical filling is that between the global minimum of 
W{£) at and the extremum at £ = oo there is no poten- 



tial barrier. Thus, walls exhibiting critical wetting nec- 
essarily form wedges exhibiting critical filling. However, 
for walls exhibiting first-order wetting, the filling tran- 
sition is first order or critical depending on whether the 
transition temperature Tp is greater or lesser than the 
spinoidal temperature T s (< T w ) at which the potential 
barrier in W(£) appears. Since the macroscopic condition 
O(Tp) = a implies that Tp can be trivially lowered by 
making the wedge angle more acute, it follows that walls 
exhibiting first-order wetting transitions will, in general, 
exhibit both types of filling transition (see Fig. 2). Note 
that the tricritical value of the wedge angle a* separat- 
ing the loci of first and second-order filling transitions 
will itself be rather small for weakly first-order wetting 
so that the Hamiltonian (Q) is still valid. The MF value 
of height critical exponent (3q for critical filling follows di- 
rectly from the equation for £q if we write the asymptotic 
decay of the binding potential as W{£) « — A£~ p where 
A is a (positive) Hamaker constant and p depends on the 
range of the forces. For systems with short-ranged forces, 
this decay is exponentially small. A simple calculation 
then yields (3q = l/p (quoted in ref . Q and implicit in ref- 
erence H) so that, for dispersion forces (corresponding 
to p = 2), the MF prediction is /3o = l/2 whilst for short- 
ranged forces /3o = 0(ln). The structure of the MF height 
profile £{x) is particularly simple near critical filling Q 
and has crucial consequences. In essence, the interface 
is flat (i.e., £(x)^£ ) for \x\<£o/a whilst for \x\>£o/a, 
the height decays exponentially quickly to its asymptotic 
planar value £ n above the wall. Importantly, the length- 
scale controlling this exponential decay is the wetting 
corr elation length £|| = y/E/W"^^) which remains mi- 
croscopic at the filling transition. One consequence of 
this is that the lateral width of the filled portion of the 
wedge is trivially identified as £ x ~2£g/a so that v x = (3$. 
More important consequences of the height structure are 
considered below. 

We now turn to the main body of our analysis con- 
cerning the nature of fluctuation effects at critical fill- 
ing and consider first fluctuations about the MF profile 
£(x) as measured by the height-height correlation func- 
tion H(x,x';y) = ( S£(x, y) 5£{x' , y') ) where 5£(x,y) = 
£(x, y) — ( £(x, y) ) and y = y'—y. To calculate the correla- 
tion function, we first exploit the translational invariance 
along the wedge and introduce the structure factor 



S(x 1 x'\ Q) = I dy e lQv H(x,x';y). 



(2) 



The assumption of MF theory is that fluctuation about 
£{x) are small and hence a Gaussian expansion of H[£] 
about the minimum suffices to determine the correla- 
tions. This leads to the differential (Ornstein-Zernikc) 
equation 

(-£<9 2 + £Q 2 +W"(£(x)-a\x\)) S(x,x';Q)=6(x-x') 

(3) 
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where we have adsorbed a factor of fc^T into the def- 
initions of £ and W(£). The structure of correlations 
across the wedge is manifest in the properties of the ze- 
roth moment Sq(x, x') = S(x, x'; 0) which can be obtained 
analytically using standard methods. We find 



S (x,x') = (\£(x)\-a) (\£(x')\-aj : 
1 H{xx') dx 



(4) 



2aW'(£ ) 



(l(x) -aX 



where H(x) denotes the Heaviside step function (H(x) = 
1 for x > 0, H(x) = otherwise). From the properties 
of the equilibrium profile £(x), it follows that the length- 
scale £ x also controls the extent of the correlations across 
the wedge. In fact, it can be seen that correlations across 
the wedge are very large and also (essentially) position 
independent, provided both particles lie within the filled 
region, implying that, at fixed y, the local height of the 
filled region fluctuates coherently. On the other hand, the 
correlations are totally negligible if one (or both) parti- 
cles lie outside the filled region since their asymptotic 
decay is controlled by the microscopic length These 
are important remarks central to the development of a 
general fluctuation theory of wedge-wetting. 

Turning next to correlations along the wedge, we note 
that a simple extension of the above analysis shows that 
the dominant singular contribution to the structure fac- 
tor has a simple Lorentzian form 



S(x, x'\ Q) 



■So(0,0) 

i + Q 2 e v 



\x\,\x'\<£ x /2, (5) 



with Sq(0, 0) —a/2W'(£o) which shows a very strong di- 
vergence as T — » Tp . The correlation length along the 
wedge is identified by £„ w (Z£ /W'(£ )) 1/2 - Substitut- 
ing for the form of W(£), and recalling the divergence 
of £o at critical filling, leads to the desired MF result 
v y = 1/p + 1/2 for the correlation length critical expo- 
nent as T — > Tp at bulk coexistence. Note that £ y ^> £ x 
so that the fluctuations are highly anisotropic and are 
totally dominated by modes parallel to the wedge direc- 
tion. The final lcngthscale that we calculate within the 
present MF/Gaussian analysis is the mid-point width £j_ 
defined by £l = {(£(0, y) - £ a ) 2 ) = H(0, 0; 0) which may 
be obtained from the Fourier inverse of S(x, x';Q). This 
leads to the intriguing relation 
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E* 



(6) 



which is one of the central results of this paper. In this 
way, we are led to the remarkable prediction that the di- 
vergence of £j_ at critical filling is universal, independent 
of the range of the intermolecular forces, and of the form 
£|| ~ i -1 / 4 which should be observable in experimental 
and computer simulation studies. We shall argue below 



that this result is not affected by fluctuation effects even 
when MF theory breaks down. 
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FIG. 2. Schematic surface phase diagram showing temper- 
ature vs. the opening angle a for a system undergoing a 
first-order wetting transition at T w in the planar case (a — 0). 
The filling transition is only first-order (Fi) if it takes place 
at a temperature above the spinoidal temperature T 3 but be- 
comes second-order (F2) if the filling temperature is less than 
T s . 

The first step in the development of a fluctuation the- 
ory for filling transitions is the derivation of a Ginzburg 
criterion. The MF analysis presented above should be 
valid if the fluctuations in the interfacial height are rel- 
atively small. Thus, we require £,±/£q < 1 or, equiva- 
lently, i 1 /^ -1 / 4 <c 1, implying that MF theory, and the 
values of critical exponents quoted above are valid in 
three-dimensions only for p < 4. For p > 4, fluctuation 
effects dominate and we can anticipate that the rough- 
ness £_l is comparable with the interfacial height £q. One 
way of approaching this problem is to formulate a renor- 
malization group theory based on the effective Hamil- 
tonian ([!]). This is an extremely difficult task and one 
which we believe is unnecessarily complicated. In view 
of the extreme anisotropy of fluctuations at filling tran- 
sitions and their coherent nature across the wedge, we 
propose that the only fluctuations that are relevant for 
the asymptotic critical behaviour are those arising from 
the thermal excitations of the mid-point height £o(y) 
along the wedge. More specifically, for a constrained 
non-planar configuration for the mid-point distribution 
{£o(y)}, we assume that all other fluctuations are small 
and, hence, following established methods jflj, may be 
treated in saddle-point approximation. Thus, we are led 
to a simpler wedge Hamiltonian (of reduced dimension- 
ality), F[£ (y)] = win' H [£(x, y)] where the dagger de- 
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notes the constraint that £(0,y) = £o(y) Vj/. In this way, 
we have derived the simpler one-dimensional model (of 
three-dimensional filling) 



F[io] = Jdy 



Q 



dy 



V F (£ ) 



(7) 



where the coefficient of the gradient term is the local 
height dependent line tension describing the bending en- 
ergy of long-wavelength fluctuations along the wedge and 
Vf is the effective wedge filling potential which has the 
general expansion 



V F {i) = 



S(6 2 - a 2 



■1+ 



A 



[p - l)a 



f-P + ■■■ 



(8) 



Note that, in the critical regime, (O(T) — a)~t, so that 
minimisation of (g) identically recovers the MF result for 
£q. For p= 1, the second term in (||) is logarithmic whilst 
for short-ranged forces, it is exponentially small. 

We propose that the effective Hamiltonian (0) contains 
all the essential physics associated with the asymptotic 
critical behaviour at filling transitions. Two checks on 
this hypothesis are that, in MF and Gaussian approxi- 
mation, the new model identically recovers the equation 
for the mid-point height and structure factor emerging 
from the more complicated model (Q) in the same ap- 
proximation. The great advantage of the new model is, 
of course, that due to its one-dimensional character, it 
can be studied exactly using transfer-matrix techniques. 
The (normalized) eigenfunctions ^n(^o) an d eigenvalues 
E n of the spectrum are found by solving the differential 
equation (setting ksT=l for convenience) 



a^(£ ) , 3a<(4>) 



E* 



(9) 



from which the quantities of interest can be calculated. 
In particular, the probability distribution for the mid- 
point height V(£q) — \ipo(£o)\ 2 and the wedge correlation 
length £ y = l/(Ei — Eo). The solution of this eigen- 
value problem for the wedge potential (^) gives a com- 
plete classification of the critical behaviour at critical fill- 
ing. The calculation confirms that MF theory is valid for 
p < 4, whilst the criticality is fluctuation dominated for 
p > 4 and is characterised by universal critical exponents 
Pa = v x = v±. = 1/4 and v y — 3/4. These exponents 
are pertinent to systems critical filling occurring in sys- 
tems with short-ranged forces and may be tested in Ising 
model simulation studies similar earlier work on critical 
wetting ||. For experimental systems with dispersion 
forces (p = 2), our predictions are (3q = v x = 1/2, v± = 1/4 
and v y = l. 

To finish our article, we make two final remarks. 
Firstly, out of bulk two-phase coexistence (5/i = fi sa t (T) — 
/i>0) and close to filling, the mid-point height, correla- 
tion lengths and roughness show scaling behaviour. For 



example, in the fluctuation-dominated regime, the solu- 
tion of (§) shows that £ = t- 1 / 4 A(<S/i*- 5 / 4 ) where A(C) is 
an appropriate scaling function. Thus, along the critical 
filling isotherm (T — Tp, dp.—>0), the height diverges as 
£o which may be easier to observe in experimen- 

tal and simulation studies. Secondly, the effective filling 
model that we have introduced can also be used to study 
complete filling occurring for T>Tp as Sfj,— >0. The crit- 
ical behaviour here is found to be MF-like (i.e. £± <C^o) 
but also universal, independent of the range of the forces 
and is consistent with the hypothesis that the geometry 
of the wedge determines the critical behaviour for this 
transition [nOl. Fluctuation effects at this transition are 
rather less interesting than for critical filling. 

In summary, we have developed a fluctuation theory 
for critical effects at three-dimensional wedge-wetting or 
filling transitions and given a complete classification of 
the possible critical behaviour. Fluctuations have been 
shown to have a much greater influence in the critical 
behaviour compared to that occurring for wetting transi- 
tions (in <i = 3) at planar surfaces and lead to a universal 
roughness exponent v±_ = 1/4. We believe that these pre- 
dictions are open to experimental verification, in wedge 
systems made form substrates exhibiting first-order wet- 
ting. 
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